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Abstract 

Consider a correlated Gaussian random energy model built by suc- 
cessively adding one particle (spin) into the system and imposing the 
positivity of the associated covariance matrix. We show that the va- 
lidity of a recently isolated condition ensuring the existence of the 
thermodynamic limit forces the covariance matrix to exhibit the Parisi 
replica symmetry breaking scheme with a convexity condition on the 
matrix elements. 
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The existence of the thermodynamic limit has been recently proved for 
the Sherrington-Kirckpatrick (SK) model fL, and more generally for any 
correlated Gaussian random energy (CGREM) model including the Derrida 
REM and the Derrida- Gardner GREM [2|. In this letter we point out that 
the proof of |2] may shed some light on the origin of Parisi's algebraic ansatz 
for the replica symmetry breaking (RSB) scheme which lies at the basis his 
solution [3] of the SK model. 

Algebraically, Parisi's ansatz may be described as follows: start from the 
one by one matrix g(0) > 0. Take q(l) as < g(0) < g(l) < 1 first replicate 
the system 
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As is well known, the RSB method consists in solving the saddle point equa- 
tions for the SK free energy after n replicas assuming that the saddle point 
matrix has the form Q , the maximum is to be found among the functions 
on [0, 1] taking the values q(k) at the prescribed points at step k of RSB. 
Performing first the limit k —>■ oo and then the limit n — > one gets: 



q(x) is continuous on [0, 1] and is the (Parisi) order parameter of the local 
square magnetization because its inverse function x(q) is such that x(q) = 



/ P(s) ds where P(q) is the overlap probability distribution between pure 
Jo 

magnetization states. Hence q(x(q)) = q and the elements q(l), ■ ■ ■ , q(k) are 
overlaps. 

Even though it has been recently proved that the Parisi free energy is 
at least a lower bound of the SK free energy [4|, the RSB is still far from 
a mathematical understanding. A major puzzle is still represented by the 
origin of the algebraic structure of the Parisi ansatz, which lies at the basis 
also of the Derrida- Gardner GREM[5|. 

In this letter we point out that the algebraic Parisi ansatz and, equiva- 
lently, the Derrida- Gardner GREM are generated by a precise prescription 
to reach monotonically the thermodynamic limit with a family of correlated 
gaussian random energy models (CGREM). 

More precisely, we will build by recurrence over N, increasing the size by 
1 at each step (equivalently, adding a spin) a family E a (N) of CGREM for 
which we require: 

1. E a (N) (iV-spins) has to be a subsystem of E a (N + 1) (N + 1-spins); 
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2. The sequence E a (N) fulfills the conditions ensuring the existence of the 
thermodynamic limit. 

The two previous conditions imply that the matrix sequence Cn fulfills the 
RSB scheme, in the sense that for each N the correlation matrix Cjy fulfills 
the Parisi ansatz at step N. In particular E a (N) can be identified with the 
Derrida-Gardner GREM with a convex growth scheme. 
For N = the system is the vacuum (no particles) with random energy 

E = £o(0) (5) 

£o is a centered Gaussian variable with variance c(0), which we can consider 
as a one-by-one matrix. We include the system of spins into a new one- 
spin system (for short: we add a spin) assuming at first that the one-spin 
Hamiltonian doesn't depend on the added spin 

£iK> = £ (0) ; (6) 

the covariance matrix of this process is 

/ c(0) c(0) \ 

c(i) = ; (7) 

V c(0) c(0) J 

This construction was called lifting in [2|. 

The one-spin system is a CGREM of size one only if its covariance ma- 
trix is a (non-degenerate) positive definite matrix. The elimination of the 
degenarcy requires the introduction in the Hamiltonian of a new Gaussian 
variable £i parametrizing the dependence on the added spin: 

£i(<n) = &(!) + 6(i)^i ; (8) 
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the covariance matrix element is 



Av(E 1 (a 1 )E 1 (r 1 )) = Av(£ (l) 2 ) + Av(^)>in ; 



(9) 



Since 



OT = 25 a)T - 1 



(10) 



we have 



Av(E 1 (a 1 )E 1 (r 1 )) = Av(^) 2 ) ~ Av(^(l) 2 ) + 2Av(Z 1 (l) 2 )6 n „ (11) 



with c(l) > c(0) 

To iterate the procedure, let us first describe the second step, i.e. the 
addition of a second spin to build a 2-spin system. As before, we first assume 
independence on the newly added spin variable. The covariance matrix turns 
then out to be (with the lexicographic order of the spin configurations) 



which we will write as 



Av(E 1 (a 1 )E 1 (T 1 )) = a (i) + a 1 (l)6 l 



(12) 



Defining c(l) = J2 k=0 <ik we have 




(13) 



C(2) 



f c(l) c(l) c(0) c(0) N 

c(l) c(l) c(0) c(0) 

c(0) c(0) c(l) c(l) 

\ c(0) c(0) c(l) c(l) , 



(14) 
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Again, this matrix is the covariance of a CGREM process only if it posi- 
tive definite and non-degenerate. As above, this requires the dependence on 
the second spin variable. Among the possible ways to parametrize this de- 
pendence we choose the minimal one, namely the preceding one which only 
modifies the subprincipal diagonals: 



(7(2) 



(15) 



( c(l) q 2 c(0) c(0) ^ 

q 2 c(l) c(0) c(0) 

c(0) c(0) c(l) g 2 

\ c(0) c(0) & c(l) y 

with c(0) < q 2 < c(l). We relabel the elements c(l) — > c(2) and g 2 — ► c(l) 
obtaining the final 2-spin covariance matrix 



C(2) 



(16) 



c(2) c(l) c(0) c(0) 
c(l) c(2) c(0) c(0) 
c(0) c(0) c(2) c(0) 

\ c (o) c (o) c (i) C (2) y 

with c(2) > c(l) > c(0) > 0. This last condition implies that the matrix is 
positive definite because the 4 principal minors are 



and 



Ai = c(2) > , 

A 2 = c(2) 2 - c(l) 2 > , 
A 3 = [(c(2) - c(l)][c(2)(c(l) + c(2)) - 2c(0) 2 ] > 

A 4 = [ C (2)-c(l)] 2 [((c(l) + c(2)) 2 -4c(0) 2 ]>0 



(17) 

(18) 
(19) 

(20) 
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Corrispondigly with the spin representation we would have 



E 2 (a 1 ,a 2 ) = £ (2) + &(2)<7i + £ 2 (2)<7 2 + 6,2^2 



(21) 



and 



Av(E 2 (a 1 , a 2 )E 2 ( Tl , r 2 )) = Av($>(2) 2 ) - Av(Ci(2) 2 ) - Av(£ 2 (2) 2 ) + Av{t lt2 {2) 2 ) 



AviE^u <7 2 )£i(7iT 2 )) = ao(2) + ai(2)(J CTliCT2 + a 2 {2)5 



In general we will assume that our construction is done adding at each step 
the iV-th spin variable and the newly added interaction terms are indipendent 
realization of the same Gaussian distribution or, in other terms, the distribu- 
tion of the Gaussian random variables depends only on the last index. The 
general scheme is the described by a correlated Gaussian process 



E N (a) =£o(N)+ &(N)<Ti+ kj(N)a l a j + ... + Ci,2,..., N a 1 ---a N 

l<i<N l<i<j<N 

(24) 



in which the distribution of each Gaussian variable depends only on the last 
index; its covariance matrix: 



+ S ntn [2Av(^(2) 2 )-2Av(^(2) 2 )] 
+ 5 CT2 , T2 [2^(6(2) 2 ) - 2Av(^ 2 (2) 2 )} 
+ <W l( W 2 4^(£ li2 (2) 2 ) ; 



and since we choose 



Av(U2) 2 ) = 2Av(^ 2 (2) 2 ) 



(22) 



the covariance matrix exibit the RSB scheme: 



N 



k 




(25) 
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turns out ot fulfill the RSB structure An important observation is that we 
may, within the positivity conditions, chose a subset of c's that fulfill the 
conditions for the existence of the thermodynamical limit as shown in [2]. 
Here we observe that those conditions translate into a convexity property 
for the matrix elements. In fact introducing the two complementary liftings 
(labelled as left I and right r) of the one-spin system: 



Q (2) 



and 



( c(2) q(2) c(0) c(0) \ 

c(2) c(2) c(0) c(0) 

c(0) c(0) c(2) c(2) 

y c(0) c(0) c(2) c(2) y 

^ c(2) g(0) c(2) c(0) ^ 



(26) 



C r (2) 



(27) 



c(0) c(2) c(0) c(2) 
c(0) c(2) c(2) c(0) 
y c(2) c(0) c(0) c(2) / 
the conditions developed in [2] are, for the symmetric sub-division of the 
system into 2 subsystems of the same size, 



C(2)<-[C,(2) + C P (2)] 



(component-wise) which imposes 



1 



c(l)<-[c(0) + c(2)] 



(28) 



(29) 



For general iV the conditions in t 2j turn out to be a set of 2^ relations between 
the matrix elements: 



Q{N)<l[Ci(N) + C r {N)] 



(30) 



With respect to the size N — 1 the only new relationships required only 
involve the quantity c(N),c(N — 1), c(N — 2). Indeed at step N the only new 
inequalities (with respect to those at step N — 1) are those coming from the 
four by four principal minors of the three matrices. After the operation of 
re-enumeration these inequalities become (element wise) 

Q(4) (A r) < i[cf } (iV) + Cf \N)} (31) 



with 



Q (4) (A0 
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c(N) c(N) 

c(N) J 



and 



/ 



C^{N) 
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\c(N-2) 

from which we have 



c(N) 
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c{N -2) 
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c(N) 
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c(N) 
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c(N -2) \ 

c(N) 
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c(N) J 



(32) 



(33) 



(34) 



c(N - 1) < -[c(N) + c(N - 2)} . 



(35) 
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Equivalently in the case of general partition of the iV-particle systems into 
two subsystems of size N\ and N2 the existence of the thermodynamic limit 
will be assured by 

c(N - 1) < ^rc(N) + ^c(N - 2) , (36) 

so that any positive increasing sequence c(i), i = 1, N with the convexity 
property 

c(i) < OLidi + 1) + (1 - aci)c(i - 1) , (37) 

for suitable chosen < ctj < 1 provides a CGREM with RSB structure whose 
thermodynamic limit is reached monotonically. 
Remark: introducing [3] the ultrametric overlap 

d N (a,r) = min(z | a { ^ r<) (38) 

one can see that the covariance matrix (|25j) is a growing function of for 
every choice of the a's. Obviously at zero temperature the overlap (|38|) co- 
incides with the standard GREM overlap among magnetization states. 
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